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Abstract 

A structure is proposed for the mass matrices of the quarks and leptons 
that arises in a natural way from the assumption that the breaking of SO (10) 
gauge symmetry is achieved by the smallest possible set of vacuum expectation 
values. This structure explains well many features of the observed spectrum 
of quarks and leptons. It reproduces the Georgi-Jarlskog mass relations and 
postdicts the charm quark mass in reasonable agreement with data. It also 
predicts a large mixing angle between and Ur, as suggested by atmospheric 
neutrino data. The mixing angles of the electron neutrino are predicted to be 
small. 



^Work supported in part by Department of Energy Grant Numbers DE-FG02-91ER-40626, DE- 
FG02-90ER-40542. 

^Email: albright.fnal.gov, "^babuOias . edu, ^smbarrSbartol . udel . edu 



In this Letter we propose a structure for the quark and lepton mass matrices that 
arises naturally in supersymmetric 5*0(10) from the simple assumption that 5*0(10) 
is broken to the Standard Model by the smallest possible set of vacuum expectation 
values (VEVs) . This structure reproduces many of the features of the known fermion 
mass spectrum. It also predicts a large value for the u^j^ — Vj- mixing angle, as is 
suggested by the atmospheric neutrino data [1]. Usually this angle is small (or not 
predicted) in grand unified theories, but in the present model its large value has a 
simple group-theoretical explanation. 

The smallest set of vacuum expectation values that can break 5*0(10) to the 
Standard Model consists of one adjoint (45) and one spinor (16) [2]. The spinor 
plays two necessary roles: it breaks the rank of the group from 5 to 4, and provides 
superlarge masses for the right-handed neutrinos. The adjoint also plays two roles: 
it completes the breaking of 5*0(10) to the Standard Model (SM) group 5*t/(3) x 
SU{2) X U{1) and produces without fine-tuning the "doublet-triplet sphtting" — 
that is, gives superlarge mass to the color-triplet partners of the SM Higgs doublets, 
while leaving those doublets light. 

Our assumption of minimality requires that there is only one adjoint Higgs. It has 
recently been shown that this is enough to break 5*0(10) with no fine-tuning, while 
preserving gauge-coupling unification [3]. Besides its economy, the postulate of having 
only one adjoint seems to be desirable in the context of perturbative heterotic string 
theory where there are limitations on multiple adjoints [4]. With only one adjoint, 
its VEV is fixed to be in the B — L direction, as this is required by the Dimopoulos- 
Wilczek mechanism for doublet-triplet splitting [5,3]. The superlarge VEV of the 
spinor is, of course, also fixed: it must point in the SU{5) singlet direction. With 
5*0(10) broken to the SM group by only these two definite VEVs, the possibihties 
for constructing realistic quark and lepton masses are quite constrained. This should 
be contrasted with other approaches that generate mass matrix textures in 5*0(10) 
utilizing extended Higgs sector [6]. 

In "minimal 5*0(10)" the quark and lepton masses come from the operators 
16jl6jl0//, where i and j are family indices, and subscript H denotes a Higgs field. 
This leads to the "naive 5*0(10) relations": N = U cc D = L, with all these matri- 
ces being symmetric. (U, D, L, and denote, respectively, the mass matrices for 
the up quarks, down quarks, and charged leptons, and the Dirac mass matrix of the 
neutrinos.) 

These relations, as is well known, lead to bad predictions: U oc D gives vanishing 
Cabibbo-Kobayashi-Maskawa angles and the relation wP^/m^ — mPg/rrii,, which is off 
by about an order of magnitude. (Superscript zero refers to parameters evaluated 
at the unification scale.) One way that U oc D can be avoided is by the quark and 
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lepton mass matrices depending on (16^^), which breaks 5*0(10). ((45//) does not 
help here, as up and down quarks have the same B — L.) However, {16h) by itself 
leaves SU{5) unbroken, which would still imply the "naive SU{5) relation" D = . 
This contains both the good prediction = m^, and the bad predictions m° = m°, 
and = m^. Therefore, the quark and lepton mass matrices must also depend 
directly or indirectly on (45//) , which is the only SU (5)-breaking VEV. Empirically, 
one finds the so-called Georgi-Jarlskog relations [7], m° = rnP^/S and m° = Sm^. Since 
(45//) (X B — L, a, natural explanation of the Georgi-Jarlskog factors of 3 and 1/3 is 
possible, as will be shown. 

The assumption of minimal VEVs for S'O(10)-breaking leads naturally, as will be 
seen, to the following forms for the quark and lepton mass matrices at the unification 
scale (with the convention that the left-handed fermions multiply them from the right, 
and the left-handed antifermions from the left): 
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These matrices, since they leave u, d, and e~ massless, are obviously not the whole 
story. At the end of this Letter, we will discuss extending the model to include the 
first generation. However, since rrie <^ m^, rrid <^ m^, and m„ <^ rric, the effects of 
such first-generation physics should be quite small on the second and third generation 
parameters that we wish to fit. It turns out that with only two parameters, e and p, 
one can get a good fit for five quantities that involve the second and third generations: 
rric/mt, nis/mb, m^/rrir, mb/rriT, and Vcb- (The other mass ratio, rrib/mt depends on 
an unknown ratio of VEVs.) 

To give some insight into the structure of the matrices of Eq. (1), and why they 
arise naturally from the assumption of minimal VEVs, it will help to explain how 
they are built up logically, layer by layer, from the heaviest generation to the lightest. 
Because the third generation is by far the heaviest, and approximately satisfies the 
SU{5) relation m° = m°, we take the first layer to come from the simple term 
I63I63IO//, giving the "1" entries in Eq. (1). 

The second-generation masses, because of the Georgi-Jarlskog factors, must de- 
pend on (45//). The simplest choice is I62I63IO//45//. This gives the "e" entries in 
Eq. (1), the factors of 1/3 just reflecting the fact that (45//) oc B — L and that a quark 
has B — L — 1/3. It can be shown that (45//) oc B — L also imphes that this term 
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contributes anti-symmetrically in flavor. (For this reason the terms 16216210^/45^^ 
and 163l63l0i/45i/ would not contribute.) 

The matrices with only the "1" and "e" entries, but without "p", would still 
not be realistic: the matrices U and D would be proportional, giving Vcb = and 
w!l/ml = m^/ml, and the Georgi-Jarlskog factor would be 9 instead of 3. (The 
"see-saw" formula would give mP^/m'^ = e^, and m°/m5 = e^/9.) 

It turns out that all three of these unrealistic features are cured in a single stroke by 
introducing a third layer that involves (16ij). The simplest term, group theoretically, 
that can be written down is of the form 162l63l6i/16'^. 16'^ is some spinor Higgs, 
distinct from 16h, which breaks the electroweak symmetry but does not participate 
in the breaking of 5*0(10) down to the Standard Model group [8]. (That is, the 
components that get VEVs arc 1(16//), and 5(16'//), where p(q) denotes a p of 
SU{5) contained in a q of 50(10).) 

This term arises most naturally from "integrating out" lO's of 5*0(10), as shown 
in Figure 1. The resulting operator is 5(162)10(163)(5(16'//))(l(16i/)). Note that 
this contributes to L and D, but not to U and N, and that it lopsidedly contributes 
to D23 and L32 but not to D32 and L23. This is the origin of the "p" entries in 
Eq. (1). This lopsidedness, which has a group-theoretical origin, explains, as will be 
seen, why the 2-3 mixing is small for the quarks {Vcb <^ 1) but large for the leptons 
(sin2 2^^, ~ 1). 

There can be a relative phase, which we will call a, between the parameters e 
and p. As is apparent from Eq. (1), this phase only enters at order e/p, which 
will presently be seen to be a small parameter. (Henceforth the symbols p and e 
will denote |p| and |e|, and the phase will appear explicitly as a.) Diagonalizing the 
matrices in Eq. (1), one finds: 

mO/mO^ l-l^(ecosa), 
K/^r = (1 - ^^(ecosa)) , 

- (1 - i^(6cos«)) , (2) 

mo/m" - eV9, 

In these expressions terms that are down by order O(e^) have been dropped. (They 
affect the results at the fraction of a percent level.) Because e is a small parameter, 
the following features of the observed masses and mixings have been reproduced by 
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the modeh the approximate equahty of m° and m°; the fact that V^°, m^^/m^, and 
m^^/m^ are all comparable, because 0(e), while ivP^/m^ is very much smaller, because 
O(e^); and the fact that m^Jm^ is about 3 times ni^/m^ (one of the the Georgi- 
Jarlskog relations). Also explained is the hierarchy among generations, which arises 
from the smallness of e and from the rank- 2 nature of the matrices. 

Since there are five observables in terms of the two parameters e and p in Eq. (2), 
the model predicts three relations among charged fermions. To study them we use 
the following input parameters: — 105.66 MeV, rrir — 1.777 GeV, ms(l GeV) = 
(180 ± 50) MeV, mb{mb) = (4.26 ± 0.11) GeV, m^{m^) = (1.27 ± 0.1) GeV [9], 
Mt = 174.1 ± 5.4 GeV and = 0.0395 ± 0.0017 [10]. The value of Mt quoted above 
corresponds to the running masses mt{mt) = 165 ± 5 GeV. 

To fit the data, various renormalization factors are needed. The factors, that will 
be denoted by rii, that run the masses from the low scales up to the supersymmetry 
scale, MsusY (taken to be at mt) are computed using 3-loop QCD and 1-loop QED 
or electroweak renormalization group equations (RGE), with inputs as{Mz) ~ 0.118, 
a{Mz) = 1/127.9 and sin^ OwiMz) = 0.2315. The relevant RGE can be found for in- 
stance in [11]. The results are (77^, 77^, 77^, 776, 77c, 77^) = (0.982, 0.984, 0.426, 0.654, 0.473, 1.0). 

The renormalization factors from MgusY up to the unification scale, Mq, are 
calculated using the 2-loop MSSM beta functions for all parameters [11], with Mq = 
2 X 10^^ GeV, and all SUSY thresholds taken to be at Msusy- These factors also 
depend on the value of tan/5, which is allowed a priori (by the perturbativity of the 
Yukawa couplings up to Mq) to be anywhere in the range 1.5 < tan/? < 65. However, 
as will be seen, within our scheme the fits constrain tan/3 to be between 10 and 
40. Results will be presented for a "central" value of 30, and where significant the 
dependence on tan/3 will be discussed. (In this model, since the light doublet, H' is a 
hnear combination of 5(10) and 5(16'), tan/3 :^mt/mb. It is also not expected to be 
very small, since the same Yukawa coupling contributes to both the top and bottom 
quark masses.) The running factors for tan/3 — 30 are {r)n/T,r)s/b,Vc/t,Vb/T,Vcb) — 
(0.956,0.840,0.691,0.514,0.873), where 77^/^- = {mym])/{mi/mj)MsusY^ and 77^6 = 

^dbl {^ch) Msusy ■ 

Aside from the running of the couplings described by the t^'s, there are finite correc- 
tions [12] to m^, mb and Vcb from gluino and chargino loops, which are proportional to 
tan/3 and thus sizable for moderate to large tan/3. These will be denoted by the factors 
(l+A^), (I+A5). and (l+Act), which depend on the supcrsymmctric spectrum: Af, ~ 

bL bp, tjj 

where f{x) = Inx/ (l—x). Ag is given by the same expression but without the chargino 
contribution (the second term) and with 6 — > s. Acb = One sees that even 
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for tan/? ^ 10, these corrections are of order 10%. The analogous corrections to 
and rrir arise only from Bino loops, while those to rric and rrit lack the tan/3 
enhancement, and so these are all negligible. 

To fit for p and e it is convenient to use the second and fifth relations of Eq. 
(2), since there is very little experimental uncertainty in m^, m,-, and Vcb- This gives 
P^[^VJim,/mr)]{^) (1-^^)(1-Ae,), and6= [^(m,/mO](^) 
(1 + e cos a p(p2^\) ) • One finds, for cosa = 1, that 

p= 1.73(1 -A,,), e = 0.136 (l-O.SAeb). (3) 

The dependence on cos a, arising only at order e/p, is rather weak: for cos a = — 1, 
p — 1.92 (1 — Acb) and e = 0.134 (1 — O.SAcb). The dependence on tan /3, because it is 
only through the renormalization factors, is also fairly weak for 10 < tan /3 < 40. For 
example, increasing tan/3 to 40 increases p by 0.7% and decreases e by 3%. Similarly, 
changing Msusy from rrit to 500 GeV only increases p by 3% and increases e by 2%. 
Henceforth, all results will be stated for tan/3 — 30, Msusy = n^t, and cos a = 1. 
Whenever results are very sensitive to these parameters, the dependence on them will 
be explicitly discussed. 

Now that p and e have been determined from Vcb and m^/m^, there are four other 
quantities that can be predicted, namely nib, and sin^ 26'j;^t-- 

(i) nib prediction: 

The first relation of Eq. (2) implies nib{nib) = m^(m^)(^^jj^^)(l — l^^^^ecoso;) 
(1 + A5). For cosa = 1, this gives nibinib) = 5.0 (1 + A5) GcV. Comparing this with 
the experimental value 4.26 ±0.11 GeV, one sees that A^ = —0.15. This is quite a 
reasonable value if tan/3 ^ 30. (With supergravity boundary conditions and a generic 
sparticle spectrum, the gluino loops contribute ~ ±0.2 to A;,, while the charginos 
contribute roughly a quarter as much and with the opposite sign [13]. We shall keep 
these numbers as a rough guide to estimate the corrections.) It should be noted that 
if tan/3 is close to 1.6 or near 60, mb{mb) will be in the acceptable range even with 
Ab = 0. However, these extreme values of tan/3 lead to wrong predictions of the 
charm mass {nic{nic) — 1.57 GeV when tan/3 ~ 1.6) and are thus disfavored within 
the model. An interesting consequence is that the model predicts the sign of p (and 
At) to be such that it decreases the 6-quark mass through the gluino and chargino 
graphs. 

(ii) nig prediction: 

The first and third relation of Eq.(2) yield ms(l GeV) = m^(m^)|e^;2^( ^ ^^^^^^ ) 
(l-iecosa^4^)(l + AJ. Forcosa = 1 this gives m,(l GeV) = 176 (1±A,) MeV. 
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Taking ~ A5 = —0.15, which is justified if the gluino contribution dominates and 
s and b are nearly degenerate, we find ms(l GeV) = 150 MeV, in excellent agreement 
with the experimental value of 180 ± 50 MeV. 

(iii) rric prediction: 

The fourth relation of Eq. (2) implies mc^m^) = mt(m()|e^(^^-^^). For cosa = 1, 
this gives rric = (1.05 ± 0.11)(1 — Act) GeV. The error refiects the la uncertainties in 
the experimental values of mt, as (= 0.118±0.004), and Vcb- (These lead, respectively, 
to 6.5%, 7%, and 4% uncertainties for mc(mc). It should also be noted that changing 
tan P by ±10 changes the rUc prediction by ^4%, changing Msusy to 500 GeV has less 
than a 2% effect, and changing cos o; to reduces rric by 3%.) Since ~ _/^j^|chargino^ 
it is reasonable to take Acb ~ —0.05, using the supergravity-like spectrum as a 
guideline. This gives rric = 1-10 ± 0.11 GeV. This is in quite reasonable agreement 
with the experimental value mc{mc) = (1.27 ± 0.1) GeV. It is interesting that the 
sign of the correction term Acb suggested by the supergravity spectrum is such that 
it improves the agreement of mc{mc) with the experimental value. 

It should also be emphasized that, whereas the predictions for and rrig were, 
in a sense, group-theoretically built into the forms given in Eq. (1), it could not be 
known in advance that the prediction for rric would come close. 

(iv) sin^ '2,6 prediction: 

The neutrino-mixing matrix Uy is defined by Vf = J2m{Ui')fm^m, where, Uf and 

are the flavor and mass eigenstates, respectively. / = e,fj,,T, and m = 1,2,3. 
Ui, — U^^^^U'^^\ where C/^^^ and U^^^ are the unitary transformations of the left- 
handed fermions required to diagonalize, respectively, L and M^, — —N'^M^^N. {Mr 
is the superheavy Majorana mass matrix of the right-handed neutrinos.) 

The crucial point, easily seen from an inspection of Eq. (1), is that a large rotation 
in the 2-3 plane will be required to diagonalize the charged-lepton mass matrix, L. 
Calling this rotation angle 62^\ one has that tan^^23'' — |-^32/-^^33| — P + ecosa. The 
actual — mixing angle is the difference between ^23'' and the corresponding 
rotation angle, ^^23''') required to diagonalize M^. 

It might appear that one can know nothing about Mi,, and therefore about 62s\ 
without knowing the precise form of Mr. However, this is not the case. Prom Eq. 
(1) it is apparent that in the limit e — > both N and M^ — —N'^M^^N are pro- 
portional to diag(0, 0, 1), so that — > 0. Thus, formally, = 0(e). If M^^ is 
parametrized by {M^%2 = {M^^)33Y/e'', and {M^^)23 = {M^%2 = {M^^^X/e, 
one finds (ignoring the first generation) that tan26'^J'^ ^ 2e\{l - X)/{1 - 2X + Y)\. 
Unless X and Y are fine-tuned, this is indeed of order e. Let k be defined by 
Re{U2s^) — Kecosa, in a phase convention where 1/23^ is real. If it is required 
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that Tnij^/iTLy^ ^ 0.05, as suggested by the atmospheric and solar neutrino data, then 
|«;| < 2. The n — t mixing angle at the unification scale is then given by 

p+ (1 - K)ecosQ; 

ia.n9nr = ^ . (4) 

1 + npe cos a 

The one-loop renormalization group equation for this quantity [14] has the simple form 
(i(lntan^^^)/(i(ln|U) = — /i^/lGvr^. Integrating yields the result that (for tan/3 = 30) 
tan6'^T- is increased by a factor of 1.03 in running down to the weak scale from the 
unification scale. 

Unlike the quark masses, the — mixing angle is very sensitive to cos a, and 
therefore sin^ 26*^^- can be in a large range, from 1 down to about 1/4. Values > 0.7 
obtain for most of the parameter range. For example, if cos a = 0, Eq. (4) simplifies to 
tan^^^ = p, giving sin^ 2^^^ = 0.78, independent of k. \i k — 0, then sin^ 2^^^ > 0.7 
for all cos a. sin^ 2^^,- reaches 1 for cos a = 1 and k — 2, and reaches ^ 1/4 for 
coscK = 1 and k — —2. (For recent attempts to generate large {v^ — Vt) niixing in 
other ways see Ref. [15].) 

There is not a unique way to extend this model to include the first generation. 
A simple possibility that gives a reasonable fit to the first-generation masses and 
mixings is to add (12) and (21) entries symmetrically to all the mass matrices. This 
would give several new predictions: (i) ^ = 3(1 + ^ecosa) (one of the Georgi- 
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Jarlskog relations), (ii) |C| = I ^Sf V^^^'l' - iV^f^I^- 

~ ^p)l- phase parameter is near tt, acceptable \Vus\ and \Vub\ 

result. The leptonic mixing angles involving the electron are given by |(f/iy)e2l — 
^(p2 + 1)1/4 + 0(e)|, |(f;^)0^| ^ I + 0(6^)1 where the 0(e) and 



0(e ) terms represent corrections from the neutrino sector. Since these mixing angles 
are both small, their precise values are sensitive to the structure of Mr. These values 
are consistent with the small angle MSW oscillations for the solar neutrinos. 

The model presented here can be tested in future experiments in several ways, (i) 
The prediction of tan/3 = 10 — 40 can be tested once supersymmetric particles are 
discovered, (ii) The spectrum of the sparticles is predicted to be such that the gluino 
and the chargino corrections to m;, decrease its value by about 15 %. (ni) More 
precise determinations of mt.,oi^{Mz) and Vcb and information about the sparticle 
spectrum will sharpen the model's prediction of 171^(111^). (iv) Large angle (z/^ — 
Vt) oscillations should be seen in long baseline experiments, but not in the ongoing 
accelerator experiments. The interpretation of the atmospheric neutrino anomaly 
in terms of (z/^ — v^) oscillations should be confirmed, (v) There are also specific 
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predictions in the model for proton decay branching ratios [16] and rare decays such 
as /i — 67 [17]. 

In this Letter we have studied a simple form for the mass matrices that is mo- 
tivated by general group-theoretical considerations, without examining a particular 
underlying unified model in great detail. That has been done in [18], where it is found 
that fermion mass matrices of the type discussed here can arise in realistic models. 
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Figure Caption 

Figure 1: A diagram that shows how vectors of fermions may be integrated out 
to produce the "p" terms in the mass matrices in Eq. (1). For group-theoretical 
reasons these produce lopsided contributions to the charged-lepton and down-quark 
mass matrices, that explain why Vcb is small while sin^ 29^r is large. 
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